Exclusive Many-Particle Diffusion in Disordered Media and 
Correlation Functions for Random Vertex Models 
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We consider systems of particles hopping stochastically on d-dimensional lattices with space- 
dependent probabilities. We map the master equation onto an evolution equation in a Fock space 
where the dynamics are given by a quantum Hamiltonian (continuous time) or a transfer matrix 
resp. (discrete time). We show that under certain conditions the time-dependent two-point density 
correlation function in the iV-particle steady state can be computed from the probability distribu- 
tion of a single particle moving in the same environment. Focussing on exclusion models where each 
lattice site can be occupied by at most one particle we discuss as an example for such a stochastic 
process a generalized Heisenberg antiferromagnet where the strength of the spin-spin coupling is 
space-dependent. In discrete time one obtains for one-dimensional systems the diagonal-to-diagonal 



transfer matrix of the two-dimensional six-vertex model with space-dependent vertex weights. For 
a random distribution of the vertex weights one obtains a version of the random barrier model 
describing diffusion of particles in disordered media. We derive exact expressions for the averaged 
two-point density correlation functions in the presence of weak, correlated disorder. 



PACS numbers: 05.40.+j, 05.60. +w, 75.10.Jm 
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I. INTRODUCTION 



Stochastically hopping particles on a lattice represent a simple model for diffusion in 
various media. Of particular interest in recent years has been the influence of static disorder 
on the dynamics. The diffusion of a single particle on a lattice with random hopping rates is 
well understood 0-0. Clearly, a system of interacting particles does not behave as a mere 
superposition of single particles. The effect of a hard core interaction on diffusion in ordered 
system has been intensively studied recently - JHJ. However the dynamics of exclusive 
particles in a random system poses a much more difficult problem. Some examples have 
been analyzed numerically [H| or analytically in a mean field approximation ||17|| . However, 
very few exact results are known. 

The main aim of this paper is to derive properties of correlation functions for certain 
diffusive systems of exclusive particles (particles with a hard-core repulsion) in a disordered 
environment. We consider d-dimensional models where the diffusion of particles is deter- 
mined by the presence of barriers of different strength. For these models the problem of 
calculating time-dependent (sometimes called "time-delayed") correlation functions in the 
steady state of the many-particle system will be shown to be reducible to a one-particle 
problem. 

We will use a description by means of stochastic occupation numbers on a lattice, the 
dynamics of which are given by a master equation in either discrete or continuous time. A 
very effective way of treating master equations is given by a Fock space method introduced 
by Doi ||18|| , further extended and applied in [I4| - f26| . In this formalism noninteracting 
particles are represented as bosons while the hard core interaction preventing the occupation 
of a lattice site by more than one particle gives rise to a fermionic character of the particles. 

Besides its usefulness for solving various problems, the mapping of the classical system 
onto a Fock space opens deeper insight into the relations to certain quantum systems. So we 
find the diffusion of exclusive particles equivalent to a Heisenberg antiferromagnet. The cor- 
relations of occupation numbers correspond to spin-spin correlations in the antiferromagnet. 
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Such a mapping is particularly interesting for one-dimensional models where the Hamilto- 
nian obtained for many equilibrium and non-equilibrium processes such as diffusion-limited 
chemical reactions or driven diffusive systems [^,0 turns out to be related to quantum 
Hamiltonians of integrable models. Thus part of the vast amount of knowledge that has 
accumulated for these models over the past decade and some of the methods applied for 
solving them such as the Bethe ansatz flH],[H are useful also for the classical systems from 
which the Hamiltonian was obtained and leads the way to new predictions. 

One obtains quantum Hamiltonians for processes in continuous time, but one may also 
study stochastic processes in discrete time. In this case the mapping to a Fock space gives 
rise to a transfer matrix which again, for many interesting problems in one dimension, turns 
out to be related to the transfer matrix of integrable two-dimensional systems such as vertex 



models P, [10| , [ll| , p7| . Here we will utilize this mapping for a study of correlated disorder in 
a one-dimensional many-particle random barrier model with exclusion. In this model the 
probability p x that a particle moves from a lattice site x to its right neighbouring lattice site 
x + 1 is the same as for a jump from x + 1 to site x. However p x is space-dependent, i.e., 
different for each pair of sites (x, x + 1). We define a parallel updating mechanism where 
these hopping rules are applied first to all pairs of sites (2x — 1, 2x) and then in a next step 
to all shifted pairs (2x, 2x + 1) as in Ref. |§. The discrete time dynamics of this model are 
thus encoded in the diagonal-to-diagonal transfer matrix of a six vertex model with certain 
randomly chosen vertex weights. 

Much is known about the (continuous-time) single-particle random barrier model with 
uncorrelated disorder, where all hopping rates p x are independently chosen random numbers 
with the same distribution g(p). The averaged probability distribution of finding the particle 
has been studied in detail by many authors 0-0. For distribution functions where the 
moments of the inverse hopping rates do not exist, a phase transition has been shown to 



take place p|fl|28[| . Only a few results are known about systems with spatially correlated 
disorder 



We will study the behaviour of the averaged density two-point correlation function in 



the steady state of the many particle system in an environment with correlated disorder. 
In Section 2 we show under which conditions this correlation function can be expressed in 
terms of the one-particle density distribution. This discussion is quite general and applies 
to ordered and disordered models in any dimension. In Section 3 we apply our result to 
systems where the corresponding quantum Hamiltonian turns out to become a generalized 
Heisenberg antiferromagnet with spatially dependent strength of the nearest-neighbour spin- 
spin coupling. We also briefly review some results on the one-dimensional case. In section 
4 we present our results on the discrete-time formulation of the process, i.e., the 6- vertex 
model. We briefly indicate the existence of a new type of phase transition in the homogeneous 
model with equal and fixed hopping rates p as they approach the value 1. Furthermore we 
derive exact expressions for the averaged correlation function for weak, correlated disorder 
and discuss the impact of the correlations on its behaviour as compared to the uncorrelated 
case. In Section 5 we summarize our results and in the Appendices we present some details 
of our calculations to Section 2 (in Appendix A) and we repeat and generalize the mapping 
of Ref. H of the diffusion problem to the 6- vertex model (in Appendix B). 

II. THE TWO-POINT CORRELATION FUNCTION FOR A MANY-PARTICLE 

SYSTEM 

In this section we will show how the two-point correlation function of a many particle 
system is related to the probability distribution of a single particle moving in the same 
environment. 

The many particle system is assumed to be described by a set of occupation numbers 
on a d-dimensional lattice n = {n,j} . Its dynamics are given by a master equation of the 
following form: 



where H' is some linear operator acting on the rij. We may also consider a discrete time 




(2.1) 
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dynamics where the time derivative in Eq. ( |2.1|) is replaced by a discrete difference. Instead 
of H' a transfer matrix T" is then used (see Section 4). 

According to Doi's formalism the master equation is mapped onto an evolution equation 
in a Fock space - 0, @: 

d t \F(t)) = H\F(t)} (2.2) 



or a similar equation in case of a the discrete time. The solution of Eq. Q2.2|) can be written 
as: 

| F(t) ) = U t \ F(0) ) with U t = e m . (2.3) 

(For a discrete time dynamics U t after t time steps is given by the i-th power T* of the 
transfer matrix T.) 

A particular configuration n corresponds to a state | n ) — J^[(c]) nj | ) where | ) is 



L 

•.t^n,■ 
3=1 

the vacuum containing no particles and L is the number of available lattice sites. A vector 



in the Fock space is decomposed as | F(t) ) = /]F(n, t) \ n ) and a scalar product is defined 

L 

by (n\m) = Yl^n^mjrijl] . 

3=1 

For noninteracting classical particles the creation operators Ca acting on site i and their 



adjoint operators Cj which annihilate particles obey bosonic commutation relations fl8|| . If 
we consider a system of particles with exclusion, i.e., if at most one particle is allowed to be 
on each lattice site, the creation and annihilation operators fulfill Pauli type commutation 
rules: operators on the same lattice site have fermionic anticommutation relations whereas 



operators for different lattice sites commute [pl] , p4]| . For simplicity we will refer to exclusive 
particles as fermions. 

Here we are interested in correlation functions. If the physical quantities A(n) and B(n) 
are analytical functions of the occupation numbers we find [pOj ]: 

< A(t)B(0) >= ( s \AU t B\ F(0) ) (2.4) 

with ( s | = ( |e^-'j Cj [18 , 2jJ and A, B are the corresponding functions of the particle 



number operators rij = c\cj. Because of the normalization condition 



-3^3- 
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(s\H = (2.5) 

which has to be valid for any H describing a stochastic process we may insert U_ t between 
the first two terms on the r.h.s. of Eq. ( |2.4| ) and obtain an expression very similar to the one 
in ordinary quantum mechanics. Note that averages are linear in the probability distribution 
here, while they are bilinear in the state in quantum theory. For special problems, however, 



linearity and bilinearity may coincide (see (|2.18|) ). 



Of particular interest are two-point correlation functions of particle numbers in a steady 
state. In what follows we restrict ourselves to systems with particle conserving dynamics, 
such as diffusive systems. A two-point correlation function for an iV-particle-system is given 
by 

G N {x, y; t) =< n x {t)n y {0) > st (2.6) 

where a subscript x denotes the label of the lattice site corresponding to a position x in space 
and n x is the particle number in x. Averaging is performed over the stationary iV-particle 
state. 

The connected two-point correlation function is defined by 

C N (x,y;t) =< n x (t)n y (0) > st - < n x > st < n y > st (2.7) 

where the arguments t and are dropped in the second term, because < n x (t) > st does not 
depend on time. 

First we consider a system with one particle only. It can be described using a master 
equation ( |2.1| ) or much easier by means of the following probability distribution: 

P(x,t;y,0) = P(x,t\y,0)P(y,0) (2.8) 

which gives the probability of finding the particle in y at time and in x at time t. The 
first expression on the r.h.s. is the corresponding conditional probability. For the two-point 
correlation function defined above we obtain 
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(1) 

G 1 (x,y;t)= Yj n x n' y F(n,t;n!_,0) 

= P(x, t\y,0)P st (y) (2.9) 

where F(n, t; rj/_, 0) is the combined probability of configuration n at time t and n at time 
0. (A superscript (k) at the sum means that the sum runs over states with a total particle 
number k only). Assuming furthermore the steady state to be homogeneous we get: 

G 1 (x,y;t) = ^P(x,t\y,0) . (2.10) 

Our aim is to express Cn in terms of P{x,t\y, 0). First we will find relations between 
Gn and G\ , one valid in the fermionic case (with particle exclusion) and another one for 
bosons. We will start with the fermionic case. The calculations for bosons are analogous 
and are not shown here for sake of brevity. 



A. The fermionic case 



It is easy to check that in the fermionic case the operators 

L L L 

V-n.,4 o 



s ,+ = E c i > S ~ = J2 4 , s> = j2(- nj + h (2.ii) 



i=i j'=i i=i 

satisfy the commutation relations of SU(2). 

In what follows we will study systems the dynamics of which are given by an SU{2)- 

symmetric Hamiltonian H, i.e., 

[#,£^ = (2.12) 
[H, S z ]- = . (2.13) 

(In the discrete time case H has to be replaced by the transfer matrix T.) 
As a consequence of this symmetry some interesting properties result: 

(i) The number of particles in the system is a conserved quantity because of Eq. ( |2.13| ). 

(ii) The Hamiltonian is hermitian. 

Furthermore, because of the normalization condition (|2.5|) and Eq. Q2.12 ) we find: 



H\N) = with \N) = —(S-) N \0) (2.14) 
and (N\H = with ( N | = — | )(S + ) N . (2.15) 

Hence, the normalized N-particle state (N < L) 

\N) norm =( L ) \N)=( L ) £ |n). (2.16) 

W W 

is a stationary solution of the problem assigning the same probability to any possible config- 
uration of the rij which has a total particle number N. (The summation in (|2.16|) runs over 

L 

states with rij = 0,1 and 52 n j = N •) The averaged occupation numbers in state | iV ) norm 

3=1 

are < n x >= N/L . 



From Eq. ( |2.6| ) and (|2.4|) we find for the correlation function: 

G f N rm (x,y;t) = ( s \n x U t n y \ N ) 

norm 

(2.17) 

L 

where the steady states ( |2.16| ) are used. Because of ( s \ = ( N |, particle conservation 

N=0 

and the orthogonality of the | ) we get: 

>,G f N erm (x,y;t) = ( N \n x U t n y \ N ) (2.18) 



Note that expression ( [2.18j ) can be understood as an imaginary time scattering matrix of a 
quantum mechanical system since the Hamiltonian is hermitian. 

Using the symmetry ( |2.12j ) the following recursion relation can be derived from Eq. ([2.18 



(see Appendix A, Eq. flA~lf)) 



, £_jV-lAf+l f IN+l 
G f ™(x,y;t) = ^—^ 7T G " m{x ^ t) + LThj ' (2 ' 19) 

From this an exact expression for GV in terms of G\ can be derived, which reduces in the 
thermodynamic limit L,N —> oo, p = N/L = const, to 

G f N erm (x, y; t) = N(l - p)G{ erm (x, y; t) + p 2 (2.20) 



where p is the mean density of the N-particle system. Inserting Eq. ( |2.1(J| ) into Eq. (2.20) 
and using the definition (|2.7|) of the connected two-point correlation function we find: 



C f N erm (x, y; t) = p(l - p)P(x, t\y, 0) . (2.21) 



Eq. ([2.21 ) is symmetric with respect to p <-> (1 — p). This reflects the particle-hole symmetry 



in a system where occupation numbers are restricted to and 1 only. The amplitude of the 
correlation function has its maximum at p = 1/2. 

Because the properties of a one-particle system are well known for many environments, 
Eq. (|2.21| ) provides a useful tool to study many-particle systems in any dimension provided 
the ,S77(2)-syinmetry holds. This will be demonstrated in the next section. 

B. The bosonic case 

Again we assume H to satisfy Eqs. fl2.12|) and (|2.13|) of the previous subsection. Note, 
however, that in the bosonic case the operators S ±,z (|2.11|) do not form a 577(2) algebra, 
instead one has a harmonic oscillator algebra with [S + , S~ ] = L where L is the number of 
sites in the system. 

As in the fermionic case the total particle number is conserved, H = is hermitian and 
Eqs. ( [2.14j ) and ( [2.15D result from the symmetry ( |2.12| ). Only the normalization constants 



and the explicit form of the steady states are different: 

AH AT! W (Jl l \ 

I iv> norm = — |AO = — E Ilr^ la>- (2-22) 

{rij} \j=i rj 

These states correspond to homogeneous probability distributions as for fermions. If the 
system is in the state | iV ) norm any configuration with particle number N has a probability 

L 

proportional to TT —. . 

In a way analogous to that in section 2.1 we find a recursion relation 



Gftg(x, y; t) = — ^— G^ son (x, y; t) + (2.23) 

from which follows 

G b r n (x, y; t) = NG b r n (x, y; t) + p 2 (2.24) 
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in the thermodynamic limit. Instead of Eq. ( |2.21| ) we obtain for the connected two-point 
correlation function for bosons: 



C°r n (x,y,z) = pP(x,t\y,0) . (2.25) 

The probability on the r.h.s. is the same as in the fermionic case because there is no 
distinction between fermion and boson for a single particle system. The connected two- 



point correlation function differs from that in the fermionic case ( |2.21| ) in the factor (1 — p) 
occuring in the latter one. This gives a quantitative description of the effect caused by the 
particle exclusion for the class of systems considered here. Of course, there is no particle-hole 
symmetry for the bosonic system. 

III. MANY-PARTICLE DIFFUSION IN A RANDOM BARRIER 

ENVIRONMENT 

Applying the results of the last section we will study many-particle diffusion in random 
environments which are constant in time. The one-dimensional system of interest is assumed 
to be described by a fermionic Hamiltonian (similar to ||14| , |2T|| ): 



3 

= o ~ 1] with °j = i c ) + c i> -i( c ] + c i)> 1 - 2n A ■ ( 3 - 2 ) 

3 

The components of <fj are the Pauli matrices. We see from Eq. ( |3.2|) that H is the Hamilto- 
nian of a generalized Heisenberg antiferromagnet with space-dependent spin-spin coupling. 
Hence the correlation functions calculated below have a physical meaning for this quantum 
system (see section 2.1). 

In what follows we will concentrate on the classical system described by the evolution 



equation ( |2.2j ) with the Hamiltonian ( |3 . 1| ) . The effect of particle exclusion is reflected in 
the terms Cj±\Cj± l and in the Pauli commutation relations between the operators. The 
Hamiltonian of nonexclusive (bosonic) particles moving in the same environment is obtained 



from (|3.1|) by omitting the aforementioned factors. 
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The Hamiltonian ( |3.1| ) generalizes the problems investigated in P,|T4jl since it allows the 
hopping rates pj to be space dependent and stochastic. The general form ( |3.1[ ) includes 
classical diffusion, diffusion in structured environments and in certain disordered media as 
e.g. in the random barrier medium. Other problems such as the random trap model are not 
described by ( ^.1[ ) and cannot be treated in the simple way proposed here. 

The Hamiltonian (|3.1| ) is SU(2) symmetric. (The corresponding bosonic Hamiltonian 
obeys the condition from section 2.2). Consequently, the relations ( [2.2 1|) and (|2.25|) between 



the correlator in the iV-particle sector and the probability distribution in the one-particle 
sector are valid in the thermodynamic limit. For a random environment ( |2.21|) and (|2.25|) are 



valid for each realization of the environment and p = N/ L is the same for each realization. 
Hence, we can average the linear equations (|2.21| ) and ( |2.25| ) over the environment (denoted 



by X). The result is: 



C f N erm (x,y;t) = p(l-p)P(x,t\y,0) (3.3) 



C b N oson (x,y;t) = pP(x,t\y,0) . (3.4) 
Eqs. Q3.3p and (|3.4j) give the averaged two-point correlation functions for the many particle 



system provided P(x, t\y, 0) is known. The latter quantity is the averaged solution of: 
d t P(x, t\y, 0) = p x [P{x + 1, t\y, 0) - P(x, t\y, 0)] + p x -x[P(x - 1, t\y, 0) - P{x, t\y, 0)] . 

(3.5) 

Eq. ( p.5|) is obviously the master equation describing the motion of a single particle in the 
environment defined by Eq. (|3.1| ). It can be formally derived as a correlation function ( p.6|) 
in the one particle sector (analogous to the procedure in section 4). 

We briefly review some known results. The most simple case is that of classical diffusion, 



i.e., with deterministic p x = D Vx, for which P(x, t\y,0) = P(x,t\y,0) is well-known from 
textbooks |TT]] . For x, t 3> 1 it is given by: 

P(x,t\y,0) = -j±=e-™ (3.6) 
11 



and its Fourier-Laplace transform S(k, oj) is 

oj + Dk z 

Note that k 2 in the denominator derives from the expansion of the Fourier transform 2(1 — 
cos k) of the lattice Laplacian to lowest order in k. Keeping also higher terms leads to an 
effective frequency dependent diffusion constant D{k) = (l — k 2 /12)D (see below). Eq. ( p.6|) 
combined with ( |3.3|) and ( |3.4|) gives the two-point correlation function for the many particle 
system in a deterministic, homogeneous environment of arbitrary dimension. The result for 
the fermionic case in one dimension is known from 



Next we consider a random barrier environment, where the hopping rates p x have the 
same probability distribution g(p) for all sites x and are not correlated. We denote the 



mean value of p by p and its variance (p — p) 2 by a 2 . (We shall neglect higher moments and 
powers of a 2 in all calculations throughout this paper.) The mean value p = D is simply 
the diffusion constant of an ordered system with g(p) = 5(p — D). 

For the single-particle problem in the case of weak, uncorrelated disorder, i.e., if all 



moments (3 m = p m exist and the p x are independent random variables, the short wavelength 



- low frequency behaviour of the Fourier-Laplace transform S(k,u) of P(x,t\y,0) is known 
- @. One introduces a generalized diffusion constant D(k,u) by writing 

1 



and finds 



S(k,u)= (3.8) 
oj + V{k, (jj)k^ 



D(k, oj) = A) + DiV^ + D 2 oj + ... 

-k 2 (E + E lv ^ +...) + ... (3.9) 



Here Dj and Ej are functions of the moments fii = (p 1 — p 1 ) L of the inverse hopping 
rates. Expanding them around D =p and neglecting higher powers of a 2 as well as higher 



moments o\ = (p — p) 1 one obtains 
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D =p- 1 » D -o 2 D- 1 (3.10) 
D X = - W 3/2 • (3.H) 

Within the framework of our approximation D2 is and Eq = Dq/12, E\ = D\/Y1 and we 
can write 

12/1 D 2X ~ 2 



D(k,u) = D ( 1 - - ) ( 1 - — (1 - -VDu) ) . (3.12) 



The prefactor 1 — k 2 /12 has its origin as above (see Eq. ( |3.7| )) in the Fourier transform of the 
lattice Laplacian and is present also in the ordered system. Up to order a 2 the uncorrelated 
disorder affects only the time dependence of the correlation function. For x, t — > 00 and 



x 2 /t fixed the averaged probability distribution P(x,t\y, 0) has the same form as ( |3.6| ) with 
a diffusion constant Dq 1 = 



In the case of strong disorder (if some of the moments (3 m do not exist) P(x, t\y, 0) cannot 
be approximated by the classical diffusion distribution [l]|| and a phase transition occurs. 



This becomes apparent in the return probability P(x, t\x, 0) J2|,|I[] which is related via Eqs. 



(3.3) and (3.4) to the autocorrelation function of the many-particle system. Its large time 



behaviour is given by P(x, t\x, 0) oc t~ a with a = 1/2 for weak disorder and a < 1/2 for 
strong disorder Jl]||. For more results on the random barrier model s. - 0. A similar 
approach can be used for studying many-particle diffusion in hierarchically structured media. 
Results for the single particle problem are given in ||2l^, |3"3| , |34T | . 

Our discussion in Sec. 2 shows that all results derived for the probability distribution of 
the single particle process are also valid for the two-point correlation functions for many- 
particle systems given by ( |3.3|) and ( |3.4|) . In the next section we are going to derive some 
new results on the effect of disorder correlations. 

IV. CORRELATION FUNCTIONS IN A RANDOM 6- VERTEX MODEL 



In this section we restrict ourselves to one space dimension only and consider discrete- 
time dynamics. It was shown by Kandel et al. that the diagonal-to-diagonal transfer matrix 

13 



of the 6-vertex model for a certain one parameter family of vertex weights describes diffu- 
sion of particles with exclusion in one dimension ||. In this mapping the time evolution 
proceeds along one diagonal of a square lattice while the space extends along the diagonal 
perpendicular to the time-diagonal (see Fig.l). They consider the standard case where the 
vertex weights do not depend on the position of the vertex in the two-dimensional lattice. 
As in the preceding section we want to study space- dependent hopping rates. In the frame- 
work of the mapping to the vertex model on a square lattice this is achieved by introducing 
vertex weights which are constant on the time-diagonal of the lattice but vary along the 
space-diagonal. For the convenience of the reader we repeat this mapping in Appendix B 
and generalize it to arbitrary, space-dependent vertex weights. 

The dynamics of the exclusion process which leads to the generalized transfer matrix of 
the 6-vertex model are defined as follows: We present the state of the system with L sites 
(L even) at time t by the quantity n(t) = {rii(t), n 2 (t), . . . ,7iz(t)} where n x (t) counts the 
number of particles on site x and can take the values or 1. The time evolution consists of 
two steps. Suppose the system is in the state n{t) with t an integer. In the first half-time 
step t — > t+1/2 we divide the chain of L sites into pairs of sites (1, 2), (3, 4), (L — 1, L). If 
both sites in a pair (2x — 1, 2x) are occupied or empty then they remain so with probability 1. 
(We exclude the possibility of particle creation or annihilation.) If there is one particle and 
one hole then the particle hops to the unoccupied site in the pair with probability p2 X -i and 
remains where it was with probability 1 — P2x-i- Note that the hopping probability in such 
a pair is the same for both directions, i.e., it does not depend on whether the particle is on 
site 2x — 1 or on site 2x. These hopping rules are applied in parallel to all pairs in the chain. 
In the second half-time step t + l/2-^t+lwe shift the pairing of the chain by one lattice 
unit such that the pairs are now (2, 3), ... (L, 1) (we assume periodic boundary conditions). 
We apply the same rules as above, but the hopping probabilities in a pair (2x, 2x + 1) are 
nowp2x- From these rules one can derive a master equation for the probability distribution 
F(n,t) Q. 

Instead of working with a master equation, we directly study the transfer matrix T which 
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encodes these hopping rules as discussed in section 2. We choose as a basis of the Fock space, 
the same basis as in section 3 where the presence of a particle corresponds to spin down and 
the absence of a particle corresponds to spin up. Then the transfer matrix is given by 

T = T even T odd = if T 2j (p 2j ) f[ T 2 ,- 1 (p 2i _ 1 ) (4.1) 

3=1 3=1 

where 

%) = l-f (4-2) 

with the Pauli matrices a x ' y,z . The time evolution operator T* for t time steps is defined as 
the kth power of T if t = k is integer and as T odd T k if t = k + 1/2. We use periodic boundary 
condition and occasionally label the spatial indices from —L/2 to L/2 — 1. (Because of the 
periodic boundary condition they are defined mod L.) The local transfer matrices Tj(pj) 
act as unit operator on all sites except on the pair (i, i + 1). Each Tj commutes with the 
generators ( |2.11 ) of SU(2), so T is symmetric under the action of SU(2) and the time- 



dependent connected two-point correlation function in the steady state 

C N (x, y;t) = (s |<ZXI N ) - P 2 ( 4 -3) 

in the sector with N = pL particles is given by the correlation function G\{x,y;t) in the 
same environment in the one particle sector as in fl2.21|) . This is true for any choice of the 
Pj and we can restrict our discussion to the one particle sector. We shall omit the index 
1 in the correlator and simply write G(x,y;t). In order to avoid boundary effects we shall 
furthermore work in the thermodynamic limit L — > oo. We denote the lattice constants in 
space and time direction by a r and a t respectively. On the square lattice one has a r = a t , 
but we shall discuss later the case a r ^ a t . The integer numbers r = R/a r and t = rj{2at) 
measure the distance R in space and r in time direction resp. in units of the respective 
lattice constants. a r and a t . Note that one full time step in the time evolution corresponds 
to t = 2a t . 

We denote by | x ) the state with the particle being on site x and ( x \ is its transposed. 
The scalar product on this space is given by ( x \ y ) = 5 X)V (S x>y is the Kronecker symbol) and 
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the unit operator is 1 = J2x I x ) ( x I • The transfer matrix A restricted to the one-particle 
sector reads 

L/2 L/2 

A = A c ™A° dd = A 2j (p 23 ) £ A 2 j-i(jP2j-i) (4.4) 

3=1 3=1 

with local transfer matrices 

Ajipj) = I j )(j I + I 3 + 1 )( 3 + 1 | -Pi(| J ) - I J + 1 ))(( j I - ( 3 + 1 I) • (4.5) 

The steady state with eigenvalue 1 of A is the vector ( a | = L~ x ^ 2 J2 X I x ) an d | a ), which 
is ( s | restricted to the 1-particle sector, is its transpose. The particle number operator n x 
is simply given by n x = | x )( x \ and the correlation function ( |4.3|) is the matrix element 
( x \ A t \ y ) = (A*)a; j3/ of A 1 . This power is defined in the same way as T l : 

\ A k if t = k 

A' = { (4.6) 

I A odd A k if t = k + 1/2 . 

With these definitions one can immediately derive a recursion relation for G(x, y; t) w.r.t. x 
and t. We have for integer values of t = k 

G(x, y; t + 1/2) = ( x |A odd A*| y ) . (4.7) 

Inserting Eq. ( |4.5| ) into this expression gives 

(1 - p x -i)G(x, y; t) + p x -iG{x - 1, y; t) x even 
(1 - p x )G(x, y; t) + p x G(x + 1, y; t) x odd. 

In the same way one obtains 

4 - p x )G(x, y; t + 1/2) + p x G(x + l,y;t+ 1/2) x even 

;i - P*-i)G(x, y; t + 1/2) + p x - X G(x -l,y,t+ 1/2) x odd . 



G(x,y;t + 1/2) 



(4.8) 



G{x,y;t + 1) 



(4.9) 



From these recursion relations which are the analogue of the recursion relation ( |3.5| ) together 
with the initial condition 
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G{x,y;0) = 5 x>y (4.10) 

one can compute the exact correlation function for arbitrary values of the local hopping 
probabilities p x . 

The simplest non-trivial case is the homogeneous model p x = const. = p. In this case the 
system is invariant under translations by two lattice units and G(x, y; t) = G(x + 2z, y + 2z; t) 
depends only on whether y is even or odd and on the distance r = x — y. For the special 
choice p = 1/2 the correlation function was computed (in a different way) by Kandel et al. 
f| (see Appendix B). 

We found that for arbitrary values of p the solution to the recursion relations (|4.8| )- ( |4.9|) 
with initial condition (|4.10| ) for x even and t integer is given by 

G(x,y;t) =p 2t 5^ 2t +'j: H") h +r ' 2 \p^\l - pf k y even (4.11) 

k=l \ k J \ k-1 J 

G(«,»; t ) = f'" ( '" 1,/2 ) f' +fr - 1,/2 )p— (1 , odd . (4.12) 

k=0 \ k J \ k J 

For x odd and t integer one finds 

G(x, y; t) = G(y, x; t) y even (4.13) 
G(x,y;t) = G(y + l,x + l;t) y odd . (4.14) 



The correlator for half-odd integer values of t is given by relations (|4.8| ). For finite densities 
p = N/L one obtains the exact connected correlation function (in the thermodynamic limit) 
by multiplying the expressions ( |4.11|) with p(l — p). Note that for p = 1/2 these expressions 



simplify considerably and we recover the result of Ref. || (see Eqs. ( |B.3|) - (B.10) in 
Appendix B). 

For large t and r (such that r 2 /t remains finite) or, equivalently, for a t = af. —>■ 0, one 



finds from (4.11) 



C N (r,t) =p(l - p)(AnDt)- 1 / 2 e- r2 / 4Dt (4.15) 
with the diffusion constant D = pj(\ — p). 
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Replacing p — > pb and taking the limit 6 — ^ 0, one recovers the Heisenberg Hamiltonian 
]2|) as H = lim^o — 1). Similarly, taking the limit t — > oo in (|4.11 ) such that bt 



remains fixed leads to the correlation function ( |3.6| ) in the continuous time formulation with 
D = p. On the other hand, for p — 1 the correlation function reduces to the 5-function 
8 r pt- Thus in this case the correlation function is invariant under the scale transformations 
r — > Ar, t — > Xt corresponding to a dynamical exponent z = 1. The crossover asp approaches 
1 (i.e. D — > oo) will be discussed in a separate publication. 

Now we study the average behaviour of the correlation function in a random environment. 
We assume all p x to be distributed in the interval < p x < 1 with the same distribution 
function g{p x ) and introduce the quantity 

A x =p x -p (4.16) 

where p denotes the mean value of the distribution g. For the sake of technical simplicity we 
choose as mean value p = 1/2 corresponding the diffusion constant D — 1. As long as p is 
not close to 1, such a choice has no qualitative influence on the averaged correlation function. 
Furthermore we shall assume that the hopping probabilities of even and odd lattice sites are 
uncorrelated, 



-iA 2 , = (4.17) 



while the correlations 



A 2x A 2y = A 2x _ 1 A 2y _ 1 = J2 h{2v)5 r , 2v (4. 18) 

V 

depend only the absolute value of the distance r = \2y — 2x\. The quantity h(0) = a 2 
appearing in the sum in the r.h.s. of Eq. ( |4.18| ) is the variance of A^. Finally, we consider 



only distributions which are sharply centered around their mean value such that higher 



moments like A 2x A 2y A 2z etc. can be neglected in a perturbative expansion of the averaged 
correlation function in its moments. 
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With the definition ( pj|) of the quantities A x we write A odd = A% dd + A odd and anal- 



ogously A even = v4g ven + A even where 7 4° dd ( even ) are the transfer matrices with all p x = 1/2. 
The averaged correlation function is averaged matrix element 



G(x, y; t) = (A% tV = { x \ (A + Ag vcn A odd + A evc Mg dd + A cvcn A odd )*| y ) . (4.19) 
Neglecting all pieces with more than two Delta matrices in this expression and using 



A even A odd = we obtain the lowest order correction to the correlation function in the 
presence of disorder 

T(x,y;t)=G(x,y;t)-G^(x,y;t) 

n n—l 



Y,J2( X |^o( A ° vcn ^o dd )^(A even Ag dd )A™- i - m | y ) + (4.20) 

1=0 m=0 



n n—l 



Y,J2( x l^o(^o venAodd )^™(^o venAodd )^r^ m l y ) • 

1=0 m=0 

Here G^(x, y\ t) denotes the correlator of the ordered system with p — 1/2 and n = t — 2. 
Some calculation shows that T(x,y;t) can be expressed in terms of a sum of three-point 
correlation functions of the ordered system: 

n 

T(x, y- 1) = 16 J2 H2v) £ ((n - m)D u (x, y- m) + D v (x, y; m)) D„(0, 0; m + \) (4.21) 

v m=0 



with 



D u (x, y;k) = (x\ (A k - A k +1 ) \y + 2u) (4.22) 
D v (x, y;k) = (x \ (A k +1/2 - A k +1 ) \ y + 2u ) . (4.23) 



Note that through the averaging, translational invariance is restored. Multiplying Eq. (|4.21|) 



by p(l — p) yields the lowest order correction to the averaged two-point correlation function. 
It is easy to compute the autocorrelation function G(0, 0; t) if h(2u) = <J 2 5 Ut o, i.e., in the 



absence of disorder correlations. From Eqs. ( J4.21 ), ( |4.22| ) and ( |B.3j ) one obtains 



G(0, 0; t) = ( 1 + 4^-^ ) G<°>(0, 0; t) . (4.24) 
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For large times this has the expected form ( |3.9|) , ( |3.10|) for D = 1 and variance 4a 2 : 



G(0, 0; t) ~ (47r J D t)- 1/2 with L> = 1 - 4a 2 . 

The behaviour of the correlation function for r / and in the presence of disorder 
correlations becomes more transparent after a Fourier-Laplace transformation. The discrete 
Fourier-Laplace transform (see Appendix B) of G^°'(x,y;t) is given by 

s ^ = T=zhzrk (425) 



while Fourier-Laplace transformation of T(x, y; t) yields 

w) = 4 1 Sm2fc( , l ~ e 2 7 ) (a 2 (l - IVT^^) - L(k, uj) S(k, u) . (4.26) 
1 — e u cos^ k K ' 



The function 



s -w/2 \ l 2 H 



L(&, w) = VI - e~ w J! M 2 ^) cos OH _ u ( 4 - 27 ) 

i/=l \ 1 + v 1 — e u J 

on the r.h.s. of Eq. (|4.26|) is the contribution of the disorder correlation function h(2u) to 
the generalized diffusion constant. Eqs. ( 4.21 ) and ( (4.26 ) are the main results of this section. 



We would like to stress that up to this point all results are exact first order contributions, 
i.e., valid for arbitrary integer values of r and t > and arbitrary values of k and u > 0. 

Now we focus on the large distance behaviour and expand H(k,u) up to order y/uj. For 
uncorrelated disorder, L(k,u) = 0, we recover the result (|3.12| ) of Refs. |5],|||| by expanding 



cos/c and exp (—u) to first order in their respective arguments. In this case the generalized 
diffusion constant does not depend on the frequency. This changes for correlated disorder. 
Assuming a decay of the form h{2v) = <r 2 exp (— |2i/|/£) one obtains for the correction 
r „ N 2 r- e-^+r 1 ) -cos2A; 

L(k, uj) « (x 2 v4J — — — = . 4.28 

V ' ; V cosh {2 y /w + 2£- 1 ) -cos 2k y ' 

For 1 <C £ <C uj^ 1 ! 2 this expression becomes 

L(k, lu) « <r 2 v^ - a 2 ^ • (4-29) 

This piece is independent of k for low frequencies k oc y/u in our expansion which considers 
only contributions of order y/uj. 
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uj 1 I 2 plays the role of a crossover length scale where the correlation function changes its 
behaviour. For 1 <C uj~ 1 / 2 <C £ one obtains a fc-dependent contribution in order y/uj: 



^,„ WV ^ = ^_. ,4.30) 

For a decay of the disorder correlations of the form h(2v) ~ a 2 \2v\~ a exp (— |2i/|/£) the 
contribution of L(k,u) to D(k,u) is small for a > 1, i.e., smaller than of order -^/a;. In 
this situation the contribution of the disorder correlations can be neglected and the system 
behaves as it was uncorrelated. 

For < a < 1 and ^y/uj the disorder contribution is larger than of order y/uj giving rise 
to a qualitative change in the frequency dependence of the diffusion constant and leading 
also to a fc-dependence. As opposed to a = with infinite correlation length the contribution 
still vanishes as uj, k 2 — > 0. 

For h(2v) = o 2 \2v\~ x (for v > 1) one finds 

L(k, uj) = ^ v^ln (1 + e" 2( " +r X) - 2e" ( " +r ^ cos 2k) (4.31) 

from which the various limiting cases can be easily derived. 

Finally we set R = ra and r = ta 2 and study the scaling limit a — > keeping R and r 
fixed. S(k, uj) becomes the well-known quantity 

S(k,u) = -±- (4.32) 
uj -\- k 

while for the first order correction ( |4.26 ) we obtain 

Ak 2 ( \ ~ 

Ti(k,oj) = — ( cr 2 — y/uj lim aL(ka, uja 2 ) ) S(k, uj) . (4.33) 

UJ -\- k \ a^O J 

From this expression we realize that the contribution of disorder correlations vanishes if their 
correlation length is finite. Furthermore, if they decay with a power law, h{2v) ~ |2z/+ l|~ a , 
then its contribution still vanishes for a > 0. Only for an exponential decay h{2v) ~ 
|2z/| _a exp (— |2z/|/£) with a = 0, and an infinite correlation length £>c<j~ 1//2 , does one obtain 
a finite contribution (we study the limit uj — ► 0). « = means that all fluctuating hopping 
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rates p2 X would be equal to some quantity p even and all P2x-i would be equal to p and 
the averaged correlation function would be an average over semihomogeneous models with 
different hopping rates at even and odd time steps. On the other hand, a = oo corresponds 
to completely uncorrelated choices of the hopping rates. However, periodic correlations of 
the type a 2 cos 2vu or a 2 J2 n ^2u,n or similar aperiodic types of behaviour also give a finite 
contribution to the averaged correlation function in the infinite time limit. 

Our results show that with decreasing correlation E(fc,a>) increases until one reaches 
a — 1. For a stronger decay only the variance a 2 of the distribution function is relevant. 



V. CONCLUSION 

We have studied systems of particles hopping stochastically on lattices of arbitrary dimen- 
sion with space-dependent hopping probabilities. For exclusion models with SU(2)-invariant 
dynamics (see ( |2.11| )- ( |2.13| )) the time-dependent two-point correlation function of occupa- 
tion numbers in the steady state is shown to be proportional to the probability distribution 
in space of a single particle moving in the same environment (see Eq. ( |2.21p ). For bosonic 
systems we obtain an analogous result (|2.25|) if the time evolution operator commutes with 



the generators of the harmonic oscillator algebra. The factor of proportionality for fermions 
is different to the one for bosons. This reflects the effect of particle exclusion, i.e., of the 
hard core repulsion. Using known results for single particle diffusion in one-dimensional, dis- 
ordered media, we obtain expressions for many-particle correlation functions for the same 
models. 

The Fock space formalism we have applied reveals relations to quantum systems. In 
particular, a classical diffusive system of exclusive particles is shown to be equivalent to a 
generalized Heisenberg antiferromagnet. 

Focussing on one- dimensional systems we study a version of the random barrier model 
with spatially correlated disorder. Its time evolution is given by the diagonal-to-diagonal 
transfer matrix of the 6-vertex model with a certain random choice of vertex weights. We 
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study the steady state and derive expressions ([4.2 1|) and ( |4.26|) for the averaged time 



dependent two-point density correlation function in the iV-particle sector in the presence 
of weak disorder. These expression are exact in the lowest order of the expansion of the 
averaged function in the moments of distribution of the hopping rates. In the language of 
the two-dimensional vertex model these are arrow-arrow correlation functions in the plane. 

We compare our results with known results for systems with uncorrelated disorder. If 
the correlation length £ of the disorder correlations is finite or if the correlations h(r) have 
infinite range but decay faster than r _1 , the contribution of the disorder correlations becomes 
negligibly small for large times and the system behaves like the system with uncorrelated 
disorder. For infinite-ranged disorder with a slower decay, h(r) ~ r~ a where < a < 1, 
the correlation function changes its behaviour even after long times and the generalized 
diffusion constant D(k,u) becomes fc-dependent even in the lowest order of the expansion. 
For uncorrelated disorder the diffusion constant depends only on u in this approximation. 
Finally we studied the infinite time limit. For a > there is no contribution from the 
correlations, but for a = 0, i.e., for a correlation that does not decay for r — > oo, we observe 
a qualitatively different behaviour. 

An interesting open question is the applicability of the ideas of Section 2 to models with 
particles of different species. 
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APPENDIX A: RECURSION RELATION FOR THE 2-POINT FUNCTION 



Starting from Eq. ( 2.18 ) for the N particle correlation function and assuming the sym- 
metry ( grp we will prove the recursion relation f gl9|) . 
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For this we will use the following equations which can be proofed easily: 

S~\ N) = (N+l)\ N + l ) (A.l) 

c t| N) = n ,\ N + l ) (A.2) 

Cj \N) = (l- nj )\N-l) (A.3) 

S+\ N) = (L-N + l)\ N- 1 ) (A.4) 

[n 3) S-]-=c] . (A.5) 

yferm 



Writing down an expression similar to fl2.18|) for and using Eq. ( |A. 1| ) it results: 



(N + l)[jG^(x,y-,t) = {N + l\n x U t n y S-\N) . (A.6) 
By means of (i) and Eq. ( |A.5|) we obtain: 

(N + l)( L )G f ™(x,y;t) = (N + l\S-n x U t n y \N) 

+ (N + l\clU t n y \N) 
+ {N +l\n x U t c\\N) . 



Because of ( |A.2|) and ( |A.3| ) and ( |A.4j ) this results in: 



N ^ J y; t) = (L - N - 1) G f ^ rm (x, y;t) + (N \U t n y \ N ) . (A.7) 

Because of the normalization condition ( |2.5| ) the second term on the r.h.s. of Eq. ( |A.7[ ) is 



equal to . Thus Eq. ( |A.7| ) results in the following recursion relation: 

G J N+l [x, y; t) = — G»(i, y; t) + ^ L _ N ■ ( A - 8 ) 



APPENDIX B: THE 6- VERTEX MODEL AS A DISORDERED DIFFUSIVE 

SYSTEM 

Here we repeat the mapping of Ref. |J of a one-dimensional diffusion problem to a 
6- vertex model and generalize it to a version of the random barrier model. Consider the 6- 
vertex model on a diagonal square lattice defined as follows: Place an up- or down-pointing 

24 



arrow on each link of the lattice and assign a non-zero Boltzmann weight to each of the 
vertices shown in figure 1. (All other configurations of arrows around an intersection of two 
lines, i.e., all other vertices, are forbidden.) The partition function is the sum of the products 
of Boltzmann weights of a lattice configuration taken over all allowed configurations. In the 
transfer matrix formalism up- and down-pointing arrows represent the state of the system 
at some given time £ (Fig. 1). Each row of a diagonal square lattice is built by M of 
these vertices. Corresponding to the M vertices there are L = 2M sites in each row. The 
configuration of arrows in the next row above (represented by the upper arrows of the same 
vertices) then corresponds to the state of the system at an intermediate time £' = £ + 1/2, and 
the configuration after a full time step t" = t + 1 corresponds to the arrangement of arrows 
two rows above. Therefore each vertex represents a local transition from the state given by 
the lower two arrows of a vertex representing the configuration on sites j and j + 1 at time 
£ to the state defined by the upper two arrows representing the configuration at sites j and 
j + 1 at time £ + 1/2. The correspondence of the vertex language to the particle picture used 
in Section 4 can be understood by considering up-pointing arrows as particles occupying the 
respective sites of the chain while down-pointing arrows represent vacant sites, i.e., holes. 

The diagonal-to-diagonal transfer matrix T acting on a chain of L sites (L even) of the 
six-vertex model with space dependent vertex weights as shown in figure 1 is then defined 

by m 

L/2 L/2 

t = n n t 2j-i = rcven T ° dd • ( ri ) 

3=1 3=1 

The matrices Tj act nontrivially on sites j and j + 1 in the chain, on all other sites they 
act as unit operator. All matrices Tj and Tj> with \j — j'\ ^ 1 commute. For an explicit 
representation of the transfer matrix we choose a spin- 1/2 tensor basis where the Pauli- 
matrix er| acting on site j of the chain is diagonal and spin down at site j represents a particle 
(up-pointing arrow) and spin up a hole (down-pointing arrow). In this basis rij — i(l — <jf) 
is the projection operator on particles on site j and sf = ±«cr|) [a x ' v ' z being the Pauli 
matrices) create (sj) and annihilate (s+) particles respectively. The matrices Tj in this basis 
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are defined by 



T 



10 

1 — pj pj 

pj 1 — pj 

1 



\ 



(B.2) 



In the particle language the matrices 3} describe the local transition probabilities of 
particles moving from site j to site j + 1 represented by the weights of the corresponding 
vertices. With these identifications the vertex model with a random choice of the numbers 
Pj in the interval < pj < 1 becomes a discrete-time version of the random barrier model. 

The transfer matrix acts in parallel first on all odd-even pairs of sites (2j — 1, 2j), then 
on all even-odd pairs. In a model with transfer matrix T = x odd T even one would start the 
time evolution at an intermediate half-odd integer time step and there will be no difference 
in the physical properties of these two systems. We assume periodic boundary conditions, 
i.e., we identify site L + 1 with site 1. 

The distinction between even and odd space coordinates becomes physically meaningful 
in the limiting case when all pj are close to, or equal to one. Suppose pj = lVj. Then 
particles which are on odd lattice sites at integer time steps move to right at a constant 
rate of two lattice units in space direction per full time step while particles move to left 
with same velocity which is the velocity of light in the system. Thus we have a system 
of non-interacting massless relativistic right- and left-movers. If the pj are not equal but 
close to 1, one expects the particles to remain relativistic but massive. The homogeneous 



massive system is studied in detail in ||32|| . This feature of the vertex model makes it more 
interesting than the continuous-time formulation by the Hamiltonian ( |3.1| ) which allows 
only for nonrelativistic diffusion. In what follows we will use the term right (left)-movers for 
particles on odd (even) lattice sites also for the nonrelativistic case. 

The time-dependent connected two-point correlation function for the homogeneous model 
with p — 1/2 was computed by Kandel et al. || for full time steps in the continuum limit 
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L — > oo. We quote their result and the corresponding expression for half-odd integer time 
intervals which are used in Section 4. One obtains with r = x — y and p = N/L: 

t integer: 



a) x, y even 

b) x odd, y even 

c) x even, y odd 

d) x, y odd 




t half-odd integer: 



a)x,yeven C N (x,y;t) = p(l - p) (]-) ( 2 * 1 j (B.O 

^ \t+(r-l)/2 y 

b) x odd, y even Cjv(a;, j/; t) = p(l - p) {^j ^ (B.8) 

c) x even, y odd Cjv(a:, y; t) = p(l - p) ( 2 * 1 j (B.9) 



1\ 2 * / 2t-l 



d)z,</odd C N (x,y;t)=p(l-p)(-) " " ■ (B.IO) 

\ZJ \t+{r-l)/2j 

Because of the distinction of right- and left-movers and full and half-time steps we have 
to define carefully Fourier and Laplace transforms of space and time-dependent correlation 



functions. We define the dynamic structure function S(k, t) as the sum of Fourier transforms 
of the two-point correlation function between right-movers and left-movers: 
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S(k,t) =Y,e lkx (C(2x,0;t) + C(2x + l,l;t)) . (B.ll) 

X 

From (p.3|) one obtains for the correlation function Cn(x, y;t) of the ordered system 
with p = 1/2 at integer times 

S (k, t) = p{\ - p) (cos k) 2t . (B.12) 

The discrete Laplace transform of a time-dependent quantity fit) is defined as the sum 

oo 

/» = Ee- wt /(*) (B.13) 
t=o 

over full time steps only. For the dynamic structure function ( p,12|) we obtain 

^■-"' l-e^co,*)' • (B ' 14) 
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1 1 p p 1 — pi— p 

FIG. 1. Allowed vertex configurations in the six-vertex model and their Boltzmann weights. 

Up-pointing arrows correspond to particles, down-pointing arrows represent vacant sites. In the 
dynamical interpretation of the model the Boltzmann weights give the transition probability of the 
state represented by the pair of arrows below the vertex to that above the vertex. 
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